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We speculate that the universe may be filled with a continuum which may be called aether. Based 
on a spherical source and sink model of electric charges, we derive Coulomb's law of interactions 
between static electric charges in vacuum by methods of hydrodynamics. A reduced form of the 
Lorentz's force law of static electric charges is derived based on a definition of electric field. 
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I. INTRODUCTION 

Coulomb's law of interactions between static electric 
charges in vacuum can be written as (refer to, for in- 
stance, 111) 
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where qi and 92 are the electric quantities of two elec- 
tric charges, r is the distance between the two electric 
charges, eg is the dielectric constant of vacuum, F21 is the 
force exerted on the electric charge with electric quantity 
92 by the electric charge with electric quantity gi, f2i 
denotes the unit vector directed outward along the line 
from the electric charge with electric quantity qi to the 
electric charge with electric quantity 92- 

The main purpose of this paper is to derive Coulomb's 
law of interactions between static electric charges in vac- 
uum by means of fluid mechanics based on spherical 
source and spherical sink model of particles. 

The motive of this paper is to seek a mechanism of 
Coulomb's law. The reasons why new mechanical in- 
terpretations of Coulomb's law are interesting may be 
summarized as follows. 

Firstly, Coulomb's law is an elementary and profound 
law in physics and play various roles in the fields of 
electromagnetism, electrodynamics, quantum mechanics, 
cosmology and thermodynamics, etc., for instance, see 
[3, 0] and the references cited there. From the point 
view of reductionism, the fundamental importance of 
Coulomb's law in all branches of physics urges the reduc- 
tionists to provide it a proper mechanical interpretation. 

Secondly, the mechanism of this action-at-a-distance 
Coulomb's law remains an unsolved problem in physics 
for more than 200 years after the law was put forth by 
Coulomb in 1785 0, % 0, H, S]. From the point view of 
reductionism, we need a satisfactory mechanical interpre- 
tation in the framework of Descartes' scientific research 
program. 

Thirdly, although the Maxwell's theory of electro- 
magnetic phenomena is a field theory [l|, the concept 
of field is different from that of continuum mechanics 



0) B B B because of the absence of a continuum. 
Thus, the Maxwell's theory can only be regarded as a 
phenomenological theory. New mechanical interpreta- 
tions of Coulomb's law may help us to establish a field 
theory of electromagnetic phenomena P, [llj . 

Fourthly, there exists some inconsistencies and inner 
difficulties in the classical electrodynamics [1, [H , l^, 3] ■ 
New theories of Coulomb's law may help to resolve such 
difficulties. 

Finally, one of the tasks of physics is the unification of 
the four fundamental interactions in the universe. New 
theories of interactions between static electric charges 
may shed some light on this puzzle. 

To conclude, it seems that new considerations on 
Coulomb's law is needed. It is worthy keeping an open 
mind with respect to all the theories of interactions be- 
tween static electric charges before the above problems 
been solved. 

Now let us briefly review the long history of mechanical 
interpretations of electromagnetic phenomena. 

Descartes was the first to bring the aether concept into 
science by suggesting that it has mechanical properties 
according to E. T. Whittakerj^. He believed that every 
physical phenomenon could be interpreted in the con- 
struction of a mechanical model of the universe. 

William Watson and Benjamin Franklin introduced the 
one-fluid theory of electricity independently in 1746 [H]- 
Henry Cavendish attempted to explain some of the prin- 
cipal phenomena of electricity by means of an elastic fluid 
in 17710. Not contented with the above mentioned one- 
fluid theory of electricity, du Fay, Robert Symmer and 
C. A. Coulomb developed a two-fluid theory of electricity 
from 1733 to 1789 [1]. John Bernoulh introduced a fluid 
aether theory of light in 17520. Euler believed that all 
electrical phenomena is caused by the same aether that 
propagates light. 

In 1821, in order to explain polarisation of light, A. J. 
Frensnel proposed an aether model which is able to trans- 
mit transverse waves. Inspired by Frensnel's luminifer- 
ous aether theory, numerous dynamical theories of elastic 
solid aether were established by Stokes, Cauchy, Green, 
MacCuUagh, Boussinesq, Riemann and William Thom- 
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son (refer to, for instance, Q). 

In 1861, in order to obtain a mechanical interpre- 
tation of electromagnetic phenomena. Maxwell estab- 
lished a mechanical model of a magneto-electric medium. 
Maxwell's magneto-electric medium is a cellular aether, 
looks like a honeycomb. In a remarkable paper published 
in 1864, Maxwell established a group of equations which 
were named after his name later. 

In a remarkable paper published in 1905, Einstein 
abandoned the concept of aether [3| • However, Einstein's 
assertion did not cease the exploration of aether, for in- 
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stance, see 

Einstein changed his attitude later and introduced his 
new concept of ether [H, I regret to admit that it 
is impossible for me to mention all the works related to 
this field in history. 

Inspired by the above mentioned works, we show that 
Coulomb's law of interactions between static electric 
charges can be derived based on a continuum mechanics 
model of vacuum and a spherical source and sink model 
of electric charges. 



II. DEFINITIONS OF SPHERICAL SOURCES 
AND SPHERICAL SINKS IN FLUIDS 

The purpose of this section is to establish a definition 
of spherical sources and spherical sinks in fluids. 

If there exists a velocity field which is continuous 
and finite at all points of the space, with the exception 
of individual isolated points, then these isolated points 
are called velocity singularities usually. Point spherical 
sources and spherical sinks are examples of singularities. 

To avoid singularities, we may introduce a mechanical 
model of spherical sources and spherical sinks with finite 
radii in fluids. 

Definition 1 Suppose there exists a hard sphere with a 
finite radius a at point Pq — {xQ,yo, zq). If the velocity 
field near the hard sphere at point P — {x, y, z) is 
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where r = ^J{x — xqY + (y — yoY ~^ i.^ ~ ^oY^ r > a, is 
the distance between the point Pq and the point P , f de- 
notes the unit vector directed outward along the line from 
the point Pq to the point P , then we call this hard sphere 
a spherical source if Q > or a spherical sink ifQ<0. 
Q is called the strength of the spherical source or the 
spherical sink. 

For convenience, we may regard a spherical sink as a 
negative spherical source. 

Suppose a static spherical source with strength Q and 
a radius a locates at the origin (0, 0, 0). In order to calcu- 
late the volume leaving the source per unit time, we may 



enclose the source with an arbitrary spherical surface S 
with radius b > a. A calculation shows that 
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where n denotes the unit vector directed outward along 
the line from the origin of the coordinates to the field 
point(a;, y, z). Equation ([3]) shows that the strength Q of 
a spherical source or spherical sink evaluates the volume 
of the fluid leaving or entering a control surface per unit 
time. 



III. FORCES ACTING ON SPHERICAL 
SOURCES AND SPHERICAL SINKS IN IDEAL 
FLUIDS 

The purpose of this section is to calculate the forces 
between spherical sources and spherical sinks in inviscid 
incompressible fluids which is called ideal fluids usually. 

Suppose the velocity field u of an ideal fluid is irrota- 



tional. then we have (refer to, for instance, [28 
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the Hamilton operator. 

It is known that the equation of mass conservation of 
an ideal fluid becomes Laplace's equation(refer to, for 
instance. 
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where (/> is velocity potential, V = + ^ 
the Laplace operator. 

Using spherical coordinates(r, 0, 1^9), a general form of 
solution of Laplace's equation Eq.® can be obtained by 
sepatation of variables asjs^l 
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where Ai and Bi are arbitrary constants, Pi{x) are Leg- 
endre's function of the flrst kind which is deflned as 



(7) 



From Eq.([2]) and Eq.®, we see that the velocity po- 
tential 0(r, 0) of a spherical source or spherical sink is a 
solution of Laplace's equation ([5|). 

The following lemma is useful in the proof of the main 
theorems in this section. 

Lemma 2 Suppose (1) the velocity field u of a fluid is 
irrotational, i.e., we have u = V</), where (p is the velocity 
potential; (2) there is an arbitrary closed surface S fixed 
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in the space without any bodies or singularities inside S; 

(3) the velocity field u is continuous in the closed surface 
S . Then, we have 

^^jj)^p<i>ndS^^^jj)^p<jyndS + jj)^pxi{u-n)dS. (8) 

where D/Dt represents the material derivative in the la- 
grangian system. 

Remark. For the proof of Lemma [21 refer to, for in- 
stance, Appendix 2 in [35| . 

Theorem 3 Suppose (l)there exists an ideal fluid (2)the 
ideal fluid is irrotational and harotropic, (3)the density 
p is homogeneous, that is^ = ^ = ^ = ^— 

(4) there are no external body forces exerted on the fluid, 

(5) the fluid is unbounded and the velocity of the fluid at 
the infinity is approaching to zero. Suppose a spherical 
source or spherical sink is stationary and is immersed in 
the ideal fluid. Then, there is a force 
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3 dt 

exerted on the spherical source or the spherical sink by 
the fluid, where p is the density of the fluid, Q is the 
strength of the spherical source or the spherical sink, a 
is the radius of the spherical source or the spherical sink, 
Uq is the velocity of the fluid at the location of the spher- 
ical source induced by all means other than the spherical 
source itself. 

Proof. Only the proof of the case of a spherical source 
is needed. Let us select the coordinates {o, x, y, z} or 
{o,Xi,X2,X3} that is attached to the static fluid at the 
infinity. 

We set the origin of the coordinates at the center of the 
spherical source. Let Si denotes the spherical surface of 
the spherical source. We surround the spherical source by 
an arbitrary spherical surface 5*2 with radius R centered 
at the center of the spherical source. The outward unit 
normal to the spherical surface Si and S2 is denoted by 
n. Let r(t) denotes the mass system of fluid enclosed in 
the volume between the surface Si and the surface iS'2 at 
time t. 

Let Fq denotes the hydrodynamic force exerted on the 
spherical source by the mass system t. Then according 
to Newton's third law, a reacting force of the force Fq 
must act on the the fluid enclosed in the mass system r. 
Let Fi = — Fq denotes this reacting force acted on the 
mass system t by the spherical source through the surface 
Si. Let F2 denotes the hydrodynamic force exerted on 
the mass system r due to the pressure distribution on 
the surface S'2. Let K denotes momentum of the mass 
system r. 

As an application of Newton's second law of motion to 
the mass system r,we have 

^^=Fi+F2 = -Fq+F2, (10) 



where D/Dt represents the material derivative in the la- 
grangian system (see, for instance, [i^, [s^, 31, [H, 33, 

In order to calculate Fq, let us calculate DK/Di and 
F2 respectively. The expressions of the momentum K 
and the force F2 are 
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where the first integral is volume integral, the second 
integral is surface integral, n denotes the unit vector di- 
rected outward along the line from the origin of the co- 
ordinates to the field point(x, y, z). 

Since the velocity field is irrotational, we have the fol- 
lowing relation 



(12) 



where (j) is the velocity potential. 

Accordin g to Ost rog radsky-Gauss theorem (see, for in- 
stance, (2i,M[3ll,[3i,[3i|) and using Eq.dlSl), we have 
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Therefore, from Eq. pTjl and Eq. (fT5)l . we have 
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Applying Lemma [2] to the second integral in Eq. (fT4| , 
we have 
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dt ]fs2 JJS2 
Putting Eq.ldni) into Eq.([Tll), we have 
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Now, let us calculate F2. Accordiiig to Lagrange- 
Cauchy integral (see, for instance, 2^ 30, 31, 33, [Ij]), 
we have 



dt 2 p 
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where f{t) is an arbitrary function of time t. Since the 
velocity u of the fluid at the infinity is approaching to 
zero, and noticing <j){t) must be of the following form 
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where Bi{t),l > are functions of time t. Thus, we have 
the following estimations at the infinity of the velocity 
field 



= O 
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where (p{x) = 0(ip{x)),x a stands for limx-,a | f{x) \ 
/il}{x) ^k,{Q<k< +00). 

Applying P7)l at the infinity and using (fT51) . we have 
I u |— > 0, d(j)/dt — > andp = Poo, where Poo is a constant. 
Thus, f{t) = Pool P- Therefore, according to (fT7|) . we 
have 



d(j) p(u ■ u) 
Using ([n|) and we have 
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Now, putting (HI]) and ^ into HU]), wc have 
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For convenience, let us introduce the following defini- 
tion 
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Since the radius R of the spherical surface S2 is arbitrary, 
we may let R to be large enough. Now, making use of 
the result (5.13) in [34|, we have 



F = 0. 

Therefore, using ([24|) . ([22]) becomes 
D 



p(j)ndS. 



(24) 



(25) 



Si 



Applying Lemma [2l (|25p becomes 



pu(u • n)dS. (26) 



For convenience, let us introduce the following defini- 
tions 

/i = ^ ^P0nd5, /2 = • n)d5. (27) 



Thus, Eg. ipS)) becomes 

Fq = /l +/2. 



(28) 



Now let us calculate the two terms in Ea. (|28|) respec- 
tively. Firstly, let us calculate the integral Ii in Ea. (p7)) . 
Since the velocity field induced by the spherical source 
with strength Q is Fq.©, then according to the superpo- 
sition principle of velocity field of ideal fluids, the velocity 
u on the surface 6*1 is 



— n + uo. 
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where n denotes the unit vector directed outward. 
Since the velocity field u is irrotational, we have 

= -^ 
47ra 

where (j>o is the velocity potential respect to Uq, i.e., Uq = 

Since the density p is homogeneous, we have 
dp 

Therefore, using Eq. (PT|) . we have 
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Suppose that dQ/dt = 0. Then, using Eg. ([30|l . we 
have 

dcj) d<j)Q 1 dQ 500 
'di~ ~dt ~ 4^~dt ~ ~dr' 
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Using Eq. ((33l) and Ostrogradsky-Gauss theorem, 
Eq.® becomes 
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We speculate that the radius a of the spherical source 
may be so small that the velocity Uq at any point of the 
spherical surface may be treated as a constant. Thus, 
Ea. ([M|) becomes 
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Now let us calculate the integral I2 in Eg. lp^ 
Noticing Eq.([29|), we have 
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For convenience, let us introduce the following defini- 
tions 



ndS, 




UgdS, 



(uo ■ n)ndS, 



(uo • n)uodS. 



Thus, Eg. ([55)1 becomes 



l2 = J1 + J2 + J3 + Ji- 
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Suppose that the radius a of the spherical source may 
be so small that the velocity Ug at any point of the spher- 
ical surface 5*1 may be treated as a constant. Therefore, 
the four integral terms in Eq. (j38p turns out to be 



(39) 
(40) 

(41) 

(42) 



Ji=0, 
J4 = 0. 

Therefore, using Eq. p9ll42)) . we have 



(43) 



at the center of the spherical source. Then Eg. ([44)) can 
be obtained following the same procedures in the proof 
of Theorem El □ 

Applying Theorem [4] to the situation that a spherical 
source or spherical sink is exposed to the velocity field of 
another spherical source or sink, we have 

Corollary 5 Suppose the presuppositions (1),(2),(3),(4) 
and (5) in Theorem\3[ are valid and a spherical source or 
spherical sink with strength Q2 is exposed to the veloc- 
ity field of another static spherical source or spherical 
sink with strength Qi, then the force F21 exerted on the 
spherical source or spherical sink with strength Q2 by the 
velocity field of the spherical source or spherical sink with 
strength Qi is 



(Qi 
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(45) 

where f 21 denotes the unit vector directed outward along 
the line from the spherical source or spherical sink with 
strength Qi to the spherical source or spherical sink with 
strength Q2, a2 is the radius of the spherical source or the 
spherical sink with strength Q2, r is the distance between 
the two bodies, V2 is the moving velocity of the spherical 
source or spherical sink with strength Q2. 

If the spherical source with strength Q2 is static in the 
aether, then Eg. (|^5)) reduces to 
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Putting Eq.(135l) and Eq.lgS]) into Eq.([2H]), we arrive at 
Eq.®. This completes the proof. □ 

Theorem[3]only considers the situation that the spheri- 
cal sources or spherical sinks are at rest in fluids. Now let 
us consider the case that the spherical sources or spheri- 
cal sinks are moving in fluid. 

Theorem 4 Suppose the presuppositions (1),(2),(3),(4) 
and (5) in Theorem\^are valid and a spherical source or 
a spherical sink is moving in the fluid with a velocity Wg, 
then there is a force 



Fq = -pg(u/ - V,) 
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3 dt 
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is exerted on the spherical source or the spherical sink 
by the fluid, where p is the density of the fluid, Q is the 
strength of the spherical source or the spherical sink, a 
is the radius of the spherical source or the spherical sink, 
Uf is the velocity of the fluid at the location of the source 
induced by all means other than the spherical source itself. 

Proof. The velocity of the fluid relative to the spher- 
ical source at the location of the spherical source is 
u/ — Vg. Let us select the coordinates that is attached to 
the spherical source and set the origin of the coordinates 



IV. A SPHERICAL SOURCE AND SPHERICAL 
SINK MODEL OF ELECTRIC CHARGES 

The purpose of this section is to establish a spherical 
source and spherical sink model of electric charges. 

It is an old idea that the universe may be filled with a 
kind of continuously distributed material which may be 
called aether (refer to, for instance, 0]). 

In order to compare fluid motion with electric fields. 
Maxwell introduced an analogy between sources or 
sinks and electric charges [sj. Recently, we derive the 
Maxwell's equations in vacuum by methods of contin- 
uum mechanics based on a continuum mechanical model 
of vacuum and a point source and sink model of electric 
charges |36|. 

In 1892 [37|, Lorentz established an electromagnetic 
theory in order to derive the Fresnel convection coeffi- 
cient. There are only two types of entities in Lorentz's 
theory: movable electrons and a stagnant aether. To 
avoid singularities, the electrons were not designed to be 
singularities as Larmor's electrons in the aether field, but 
were extremely small hard spheres with a finite radius. 

Inspired by Maxwell and Lorentz 37|, we speculate 
that electric charges may not be singularities, but may be 
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extremely small hard spherical sources or spherical sinks 
with finite radii. 

In [36], we introduced a hypothesis that the constitu- 
tive relation of the aether satisfies 



dt 



1 

2^ 



1 dsjj 
^ 2G dt 



(47) 



where s^- is the stress deviator, is the strain deviator, 
G is the shear modulus, rj is the dynamic viscocity. 

Now let us introduce the following definition of 
Maxwelllian relaxation time r 

V 

Therefore, using Ea. (|48l) . Eq. (|47)) becomes 



(48) 



T 



ds 



dt 



de 
dt 



(49) 



Let Tq be the characteristic time of a observer of an 
electric charge. We may suppose that the observer's time 
scale To is very large compares to the the Maxwelllian re- 
laxation time r. So the Maxwelllian relaxation time r is 
a relatively a small number and the stress deviator Sij 
changes very slowly. Thus, the second term in the left of 
Ea. p5)l may be neglected. Therefore, the observer con- 
cludes that the aether behaves like the Newtonian-fluid. 
According to this observer, the constitutive relation of 
the aether may be written as 

s,, = 2^^. (50) 

We see that Eg. ((50)) is the constitutive relation of a 
Newtonian-fluid. Therefore, the observer concludes that 
the aether behaves as a Newtonian-fluid under his time 
scale. 

Now let us introduce the following hypothesis. 

Hypothesis 6 Suppose all the electric charges in the 
universe are small hard spherical sources or spherical 
sinks with finite radii in the aether. We define a spherical 
source as a negative electric charge. We define a spheri- 
cal sink as a positive electric charge. The electric charge 
quantity of an electric charge is defined as 

Qe = -kqpQ, (51) 

where p is the density of the aether, kq is a positive 
dimensionless constant, Q is called the strength of the 
spherical source or spherical sink. 

A calculation shows that the mass m of a electric 
charge is changing with time as 



dm Qe 
— = -pQ = — , 
at kg 



(52) 



where Qe is the electric charge quantity of the electric 
charge. 

We may introduce a hypothesis that the masses of elec- 
tric charges are changing so slowly relative to the time 
scaler of human beings that they can be treated as con- 
stants approximately. 



V. DERIVATION OF COULOMB'S LAW OF 
INTERACTIONS BETWEEN STATIC ELECTRIC 
CHARGES IN VACUUM 

The purpose of this section is to derive Coulomb's law 
of interactions between static electric charges in vacuum. 

It is known that the motion of an incompressible 
Newtonian-fluid is governed by the Navier-Stokes equa- 
tions (refer to, for instance, la, S El, [13, Hi), 



du 

'dt 



-l-(u- V)u 



-Vp - ryV^u, 



(53) 



where u is the velocity field of the fluid, p is the pressure 
field, p is the density field, rj is the dynamic viscocity 
coefficient, t is time. 

The definition of Reynolds number Re of a fluid field 

is 



Re 



(54) 



where po is the characteristic density, Uq is the charac- 
teristic velocity, Lq is the characteristic length, r/o is the 
characteristic dynamic viscocity coefficient. 

Hypothesis 7 We speculate that the characteristic ve- 
locity Uq of an electric charge is so high compares to 
the characteristic dynamic viscocity coefficient rjo of the 
aether that the Reynolds number Re of the fluid field is 
a large number. 

Under this hypothesis, we may treat the aether as an 
inviscid incompressible fluid when we study the motion 
of electric charges. Therefore, according to Hypothesis [3 
the motion of the aether is gov erned by the Euler equa- 
tions (refer to, for instance. |29l. 30l . 31, H, [13, 0] ) , 



du 

'at 



(u- V)u 



-Vp. 



(55) 



Now based on Hypothesis [5] and Hypothesis [71 we can 
make use of Theorem |3] and Theorem 2] to study the mo- 
tions of electric charges. 

Theorem 8 Suppose a static electric charge with an 
electric charge quantity q2 is exposed to the electric field 
of another static electric charge with an electric charge 
quantity qi, then the force F21 exerted on the electric 
charge with electric charge quantity 92 by the electric field 
of the electric charge with electric charge quantity qi is 



F21 = 



1 



91 92 



Skskg 47reo 47rr^ 
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TC21 



Anal d 
ikq dt 
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r2i, 



(56) 



where f 21 denotes the unit vector directed outward along 
the line from the electric charge with electric charge quan- 
tity qi to the electric charge with electric charge quantity 
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(72, 0,2 is the radius of the electric charge with electric 
charge quantity q2, r is the distance between the two elec- 
tric charges, kg and ks are two positive dimensionless 
constants. 



Proof. From Hypothesis [SI we have 
Qi = --r — I Q2 



02 



7^, (57) 



where Qi and Q2 are the strengthes of the electric charges 
respectively. From a definition in Eq.(65) in [36| . we have 



eo = 



kfi ' 



(58) 



where eg is the dielectric constant of vacuum. Putting 
Eq.dnZl) and Eq.dSS]) into Eq.®, we get Eq.dSH). □ 
If we ignore the second term in Eq. (|56p . then we get 
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(59) 



Compare Eq. ((59| with Eq. ([T]) , it is natural to introduce 
the following hypothesis. 

Hypothesis 9 Suppose we have the following relation 



kskg = -. 



(60) 



Now based on Hypothesis [H from Eq. ([59]) we see that 
we have arrived at Coulomb's law ([T]) of interactions be- 
tween static electric charges in vacuum. 

Theorem [5] only states the force exerted on a static 
electric charge by the electric field of another static elec- 
tric charge. We may generalize this result to the case of 
an static electric charge exposed to any electric field of 
the aether. 

In [3^, we introduced the following definition 



E : 



-fc, 



du 

' 9t ' 



(61) 



where u is the displacement of the visco-elastic aether, 
^ is the velocity field of the aether, E is the electric 
field intensity, fc^; is a positive dimensionless constant. 

Since the observer of an electric charge concludes that 
the aether behaves as a Newtonian-fluid under his time 
scale, we may define the electric field intensity as the 
velocity field of the aether. Therefore, we may introduce 
the following definition. 

Definition 10 We define the electric field intensity in 
the aether as 



E = -ksVa, 



(62) 



where Va is the velocity field of the aether at the location 
of an testing electric charge induced by all means other 
than the testing electric charge itself E is the electric 
field intensity, ks is a positive dimensionless constant. 



Theorem 11 Suppose a static electric charge with an 
electric charge quantity q is exposed to an electric field 
E of the aether, then the force F exerted on the electric 
charge by the electric field E is 



qF, 



4na^eo 5E 
Sfcg dt ' 



(63) 



where q is the electric charge quantity of the electric 
charge, a is the radius of the electric charge, eg is the 
dielectric constant of vacuum, kg is a positive dimen- 
sionless constant, E is the electric field at the location of 
the electric charge induced by all means other than the 
electric charge itself. 



Proof. From Hypothesis [6l we have 

kqp 



(64) 



where Q is the strength of the electric charge. Putting 
Eg. lp^ and Ea. ([M)) into Eq.® and making use of 
Eq.® and Eq.dSni), we get Eq.dM]). □ 

If we ignore the second term in Eq. (j63p , then we get 

F = gE. (65) 

We know that the Lorentz's force law can be written 



as 



F = (jE + X B, 



(66) 



where q is the electric charge quantity of a electric charge, 
F is the force exerted on the electric charge, Ve is the 
velocity of the electric charge, E and B are the electric 
field intensity and the magnetic induction respectively at 
the location of the electric charge induced by all means 
other than the electric charge itself. 

We see that Eq. (p5)) is a reduced form of the Lorentz's 
force law ([55]) for static electric charges in vacuum. 



VI. SUPERPOSITION PRINCIPLE OF 
ELECTRIC FIELD OF STATIC ELECTRIC 
CHARGES 

From the definition [TOl we see that the electric field 
intensity is a linear function of the velocity field of the 
aether. Therefore, the superposition principle of electric 
fields of static electric charges in vacuum is deduced from 
the superposition theorem of the velocity field of fluids. 



VII. CONCLUSION 

We show that the forces between spherical sources and 
spherical sinks in ideal fluids coincide with Coulomb's 
inverse-square-law of interactions between static electric 
charges. It is an old idea that the universe may be filled 
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with a kind of continuously distributed material which 
may be called aether. Based on a continuum mechanical 
model of vacuum and a spherical source and spherical 
sink model of electric charges, we derive Coulomb's law 
of interactions between static electric charges in vacuum. 
Furthermore, we define the electric field intensity as a lin- 
ear function of velocity field of the aether at the location 
of an testing electric charge induced by all means other 
than the testing electric charge itself. Thus, we obtain a 
reduced form of the Lorentz's force law for static electric 
charges in vacuum. 



VIII. DISCUSSION 

There exists some interesting theoretical, experimental 
and applied problems in the fields of fluid mechanics, con- 
tinuum mechanics, classical electrodynamics, quantum 
electrodynamics and other related fields involving this 
theory of Coulomb's law of interactions between static 
electric charges. It is an interesting task to generalize 
this work further to describe the forces exerted on a elec- 
tric charge moving with a velocity in a electromagnetic 
field. 
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